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We investigate the effect of a variation of fundamental constants on primordial element production 
in Big Bang nucleosynthesis (BBN). We focus on the effect of a possible change in the nucleon- 
nucleon interaction on nuclear reaction rates involving the A ^ 5 (^Li and ^He) and A = 8 (^Be) 
unstable nuclei. The reaction rates for ^He(d,p)^He and ^H(d,n)^He are dominated by the proper- 
ties of broad analog resonances in ^Li and ^He compound nuclei respectively. While the triple- 
alpha process '^}ie{aa,y)^^C is normally not effective in BBN, its rate is very sensitive to the 
position of the "Hoyle state" and could in principle be drastically affected if '^Be were stable dur- 
ing BBN. We found that the effect of the variation of constants on the ^He(d,p)^He, ^H(d,n)'*He 
and '*He(aa, 7) '^C reaction rates is not sufficient to induce a significant effect on BBN, even with 
a stable ^Be. The main influences come from the weak rates and the A =2, n(p,7)d, bottleneck 
reaction. 
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1. Introduction 

Constraints on the possible variation of fundamental constants are an efficient method of test- 
ing the equivalence principle [|]], which underpins metric theories of gravity and in particular gen- 
eral relativity. These constraints are derived from a wide variety of physical systems and span a 
large time interval back to Big Bang Nucleosynthesis (BBN). Using inputs from WMAP for the 
baryon density BBN yields excellent agreement between the theoretical predictions and astro- 
physical determinations for the abundances of D and ^He despite the discrepancy between the 
theoretical prediction of ^Li and its determined abundance in halo stars. The effects of the variation 
of fundamental constants on BBN predictions is difficult to model. However, one can proceed in 
a two step approach: first by determining the dependencies of the light element abundances on the 
nuclear parameters and then by relating those parameters to the fundamental constants, following 
our earlier work [||]. 

It is well known that, in principle, the mass gaps at A =5 and A = 8, prevent the nucleosynthetic 
chain from extending beyond ^He. The presence of these gaps is caused by the instability of ^He, 
^Li and ^Be which are respectively unbound by 0.798, 1.69 and 0.092 MeV with respect to neutron, 
proton and a particle emission. Variations of constants will affect the energy levels of the ^He, ^Li, 
^Be and '^C nuclei [^, 0], and hence, the resonance energies whose contributions dominate the 
reaction rates. In addition, since ^Be is only slightly unbound, one can expect that for even a small 
change in the nuclear potential, it could become bound and may thus severely impact the results of 
Standard BBN (SBBN). It has been suspected that stable ^Be would trigger the production of heavy 
elements in BBN, in particular that there would be significant leakage of the nucleosynthetic chain 
into carbon [^]. Indeed, as we have seen previously |^, changes in the nuclear potential strongly 
affects the triple-alpha process and as a result, strongly affect the nuclear abundances in stars. 



2. Thermonuclear reaction rate variations 

It would be desirable to know the dependence of each of the main SBBN reaction rates to 
fundamental quantities. This was achieved in Ref. [|5|], but only for the first two BBN reactions: the 
n-f->p weak interaction and the p(n,7)d bottleneck. Here, we propose to extend this analysis to the 
^H(d,n)'*He and ^He(d,p)^He reactions that proceed through the A = 5 compound nuclei ^He and 
^Li, and to the ^He(aa, yY^C reaction that could bridge the A = 8 gap. 

The weak rates that exchange protons with neutrons can be calculated theoretically and their 
dependence on Gf (the Fermi constant), Qnp (the neutron-proton mass difference) and (the 
electron mass) is explicit [^. The dependence of the n-i-p— ^d-i-yrate [ [Io| ] cannot be directly related 
to a few fundamental quantities as for the weak rates, but modeling of its dependence on the binding 
energy of the deuteron Bd has been proposed [11 , |I2]. 



For the ^H(d,n)'^He, ^He(d,p)'^He and '*He(o;o;, 7)'^C reactions, we used a different approach. 
In these three reactions, the rates are dominated by the contribution of resonances whose proper- 
ties can be calculated within a microscopic cluster model. The nucleon-nucleon interaction V(r) 
depends on the relative coordinate and is written as: 

V(r)=Vc(r) + (l + 5^JVw(r), (2.1) 
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where Vc(r) is the Coulomb force and V;v(r) the nuclear interaction. The parameter 5^^ charac 



terizes the change in the nucleon-nucleon interaction. When using the Minnesota force []13|], it is 
related to the binding energy of deuterium by ABo/Bd = 5.7701 x 5^^^ (The variation of the 
Coulomb interaction is assumed to be negligible compai^ed to the nuclear interaction). The next 
important step is to relate ABj) to the more fundamental parameters. To summarize, Bd has been 
related, within an co and a mesons exchange potential to quark masses and Aqcd by Flambaum & 
Shuryak [ 14 ] and subsequently to more fundamental parameters (see Coc et al. |^ and references 



therein), and in particular to the fine structure constant. 

2.1 The triple-alpha 

The triple-alpha reaction is a two step process in which, first, two alpha-particles fuse into 
the ^Be ground state, so that an equilibrium {2a -H-^Be) is achieved. The second step is an- 
other alpha capture to the Hoyle state in ^^C. In our cluster approximation the wave functions 
of the ^Be and '^C nuclei are approximated by two and three-cluster wave functions involving 
the alpha particle, considered as a cluster of 4 nucleons. It allows the calculation of the vari- 
ation of the ^Be ground state and '^C Hoyle state w.r.t. the nucleon-nucleon interaction, i.e. 
d^. In Ref. [|], we obtained £g.^.(^Be) = (0.09208 - 12.208 x 5^^) MeV, for the **Be g.s. and 
£■^(12^) = (0.2877 - 20.412 x 5^^) MeV, for the Hoyle state. From these relations, it is possible 
to calculate the partial widths, and subsequently the '*He(aa, 7)'^C rate as a function of d^^ [Q]. 
Indeed, variations of 5^^ of the order of 1%, induces orders of magnitude variations of the rate 
[^] at temperatures of a few 100 MK. In addition, one sees that ^ (^Be) (relative to the 2-a 
threshold) becomes negative (i.e. ^Be becomes stable) for 5^^^ ^ 7.52 x 10^^. In that case, we have 
to calculate the two reaction rates, ^He(a,7)^Be and ^Be(a,7)^^C for a stable *^Be. The calcula- 
tion of the rate of the second reaction can be achieved using the sharp resonance formula with the 
varying parameters of the Hoyle state from Ref. [^. For the first reaction, '*He(oc, 7)^Be, we have 
performed a detailed calculation following [ ]T5| ] to obtain the astrophysical ^-factor, and reaction 
rate, for values of the ^Be binding energy of B^ = — S^.v^^Be) = 10, 50 and 100 keV. 

2.2 The 3He(d,p)^He and 3H(d,n)4He reactions 

The ^He(d,p)^He and ^H(d,n)^He reactions proceeds through the ^Li and ^He compound nu- 
clei and their rates are dominated by contributions of | analog resonances. The corresponding 
levels are well approximated by cluster structures (^He(8)d or t(X)d), so that we can use the same 
microscopic model as for the ^He(aa,7)^^C reaction. However, unlike in the case of ^Be, the 
^He and ^Li nuclei are unbound by -^1 MeV and the resonances are broad. Therefore the issue of 
producing A = 5 bound states, or even a two step process, like the triple-alpha reaction is irrelevant. 

To be consistent with our previous work, we want to reproduce, for 5^,^=0, the experimental 5- 
factors (see references in Ref. [|l6|]) obtained by a full /^-matrix analysis, but, here, for convenience, 
we restrict ourselves here to the single pole /^-matrix approximation which will be shown to be 
sufficient: 

cT(F) (^^)' r.n(£)ro„t(£) 

^ (£*+A£*-£)2 + r2(£)/4 ^ ■ ' 

For the ^H(d,n)'*He reaction, we use the parameterization of Barker [^, which reproduces the 
resonance corresponding to the |^ state at 16.84 MeV which is in perfect agreement with the full 
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Figure 1: 5-factor curves for the ^He(d,p)^He reaction from Ref. [|l6| (dashed), from our fit (soHd and 
overlapping) and for extreme variations of 5^^ = ±0.15 (dotted). Deviations with experimental data at very 
low energy are due to screening. 



/?-matrix fit that we used in previous work [Q |8|, |T9|, 20 1. For the ^He(d,p)^He reaction, we 
performed a fit (see Fig. |T]) of the full R-matrix ^-factor provided from Ref. [16] since, in that case, 
the parameterization of Barker [17], does not reproduce well the more recent data. 

The additional level shifts obtained with our cluster model are given by AEr = —0.327 x 8^^ 
for 3H(d,n)4He and AEr = -0.453 x 5^^ for 3He(d,p)4He (units are MeV) These energy de- 
pendences are much weaker (~ 20-30 keV for \5^^ \ < 0.03) than for *^Be and ^^C. This is expected 
for broad resonances which are weakly sensitive to the nuclear interaction. In contrast, Berengut 
et al. 1^ find a stronger energy dependence. These authors perform Variational Monte Carlo cal- 
culations with realistic N-N interactions, which provide better D and ^H/^He wave functions, but 
which are not well adapted to broad resonances, such as those obsei^ved in ^He and ^Li. Besides, 
the calculation of AEr as a function of 5^^ is obtained by the difference between the energy of the 
J States and the thi^esholds for the two-clusters emission in the entrance channel, both depend- 
ing on the N-N-interaction. Berengut et al. [||] assume that these levels follow the dependence 
of the ^Li and ^He ground states, but the | resonant levels state have indeed a ^He^d or t®A 
structure, different from the ground states. We also use a more elaborate parameterization of the 
cross-section. 



3. Effects on primordial nucleosynthesis 

The results of the former sections can be implemented in a BBN code in order to compute the 
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primordial abundances of the light elements as a function of 8^^. For ^He, D, ^He and ^Li, we 
found that the effect of the ^He(d,p)^He and ^H(d,n)^He rate variations was negligible compared to 
the effect of the nop and n(p,7)d reaction rate variations that we considered in our previous work 
[^]. Hence, next, we allow those two last reaction rates to vary through the coupled variation of 
5nn' ^Dj electron and quark masses, Gp, Qnp, Aqcd> etc.... as done in Ref. Then, with updated 
D and ^He primordial abundances abundances deduced from observations, we obtained [ pT| ] 

-0.0025 < < 0.0006. (3.1) 

for typical values of the parameters. Those allowed variations in are too small to reconcile ^Li 
abundances with observations, where « —0.01 is required. We can easily extend our analysis by 
allowing both Tjio and 5^^^ to vary. This allows one to set a joint constraint on the two parameters 5^^^^ 
and baryonic density, as depicted on Figure ^ No combination of values allow for the simultaneous 
fulfilment of the "^He, D and 'Li observational constraints. 

Note that the most influential reaction on 'Li is surprisingly [||, |^ n(p,7)d as it affects the 
neutron abundance and the 'Be destruction by neutron capture. The dependence of this rate to Sd 
that we used comes from Dmitriev et al. [^] but, very recently, this has been challenged by the 



work of Carrillo et al. [12| that provide a very different dependence. If so, the influence of 5^^ on 
'Li would have to be re-evaluated. 

Finally, we investigated the production of ^^C by the ^He(aa, 7)^^C, or the ^He(a,7)^Be 
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and ^Be(a,7)^^C reactions as a function of 8,^^. This is to be compai^ed with the CNO (mostly 
^^C) SBBN production that has been calculated in a previous work [ |20| ] to be CNO/H = (0.5 — 
3.) X 10^^^, in number of atoms relative to hydrogen. A network of f«400 reactions was used, 
but the main nuclear path to CNO was found to proceed from ^Li(n,7)^Li(a,n)^^B, followed by 
"B(p,7)^2(-^ ^'B(d,n)^2(-^ iiB(d,p)i2B and ^^B(n,7)^2g reactions. To disentangle the ^^C produc- 
tion through the "^He— 7-^Be^^^^C hnk, from the standard ''Li^^^Li^'^B— ^'^C paths, we reduced 
the network to the reactions involved in A < 8 plus the '*He(o:a, 7)'^C, or the ^He(a,7)^Be and 
^Be(a,7)^^C reactions, depending whether or no ^Be would be stable for a peculiar value of 5^^^^. 
The carbon abundance shows a maximum at 5^^^ 0.006, C/Hss 10^^' [pi]], which is six orders 
of magnitude below the carbon abundance in SBBN [pO|]. This can be understood as the baryon 
density during BBN remains in the range 10^^ to 0.1 g/cm^ between 1.0 and 0.1 GK, substantially 
lower than in stars (e.g. 30 to 3000 g/cm^ in stars considered by Ekstrom et al. [Q). This makes 
three-body reactions Uke ^He(o:a,7)'^C much less efficient compared to two-body reactions. In 
addition, while stars can produce CNO at O.IGK over billions of years, in BBN the optimal tem- 
perature range for producing CNO is passed through in a matter of minutes. Finally, in stai^s, 
'^He(aa, 7)^^C operates during the helium burning phase without significant sources of ^Li, d, p 
and n to allow the ^Li— >^Li— )-^'B— )'^^C, A=8, bypass process. 

Note that the maximum is achieved for 5^^ 0.006 when ^Be is still unbound so that contrary 
to a common belief, a stable ^Be would not have allowed the buildup of heavy elements during 
BBN. This is illustrated in Figure ^ which displays the evolution of the ^^C and ^Be mass fractions 
as a function of time when ^Be is supposed to be bound by 10, 50 and 100 keV (solid lines). They 
both increase with time until equilibrium between two a-particle fusion and ^Be photodissocia- 
tion prevails as shown by the dotted lines. For the highest values of B^, the ^Be mass fraction 
increases until, due to the expansion, equilibrium drops out, as shown by the late time behavior 
of the upper curve (Bg = 100 keV) in Figure ||. For ^ 10 keV, the '^C production falls well 
below, out of the frame, because the ^Be(a,7)^^C reaction rate decreases dramatically due to the 
downward shift of the Hoyle state. For comparison, the SBBN sa400 reactions network (essentially 
the ^Li^^^Li— >"B^>'^C chain) result [20] is plotted (dashed line) in Figure^ It shows that for the 
^He— )-^Be^^^C path to give a significant contribution, not only ^Be should have been bound by 
much more than 100 keV, but also the ^Be(a, 7)'^C rate should have been much higher in order to 
transform most ^Be in '^C. 



4. Discussion 

We have investigated the influence of the variation of the fundamental constants on the pre- 
dictions of BBN and extended our previous analysis [||]. Through our detailed modeling of the 
cross-sections we have shown that, although the variation of the nucleon-nucleon potential can 
greatly affect the triple-a process, its effect on BBN and the production of heavier elements such 
as CNO is typically 6 orders of magnitude smaller than standard model abundances. Even when 
including the possibility that *^Be can be bound, at the temperatures, densities and timescales as- 
sociated with BBN, the changes in the ^He(aoi, 7)'^C and ^Be(a,7)^^C reaction rates are not 
sufficient. We have also extended our previous analysis by including effects involving ^He and 
^Li. This allowed us to revisit the constraints obtained in Ref. iQ] and in particular to show that the 
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'Be 10, 50 and 100 keV bound 




Time (s) 

Figure 3: '^C and ^Be mass fractions as a function of time, assuming ^Be is bound by 100, 50 and 10 keV 
as shown by the upper to lower solid curves respectively [^. (Only the '^C mass fraction curve, for Bg = 
10 keV, is shown; others are far below the scale shown). The dotted lines correspond to the computation at 
thermal equilibrium and the dashed line to the SBBN [ ^0| | production. 



effect of the ''He(d,p)'^He and -''H(d,n)'^He cross-sections variations remain small compared to the 
n(p,7)d induced variation. 
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